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Abstract 

The  line  mass  and  implosion  energy  inferred  from  the 
PRS  implosion  time  and  load  current  depend  on  the 
shapes  of  gas-puff  radial-density  and  current-history  pro¬ 
files.  Snowplow  calculations  with  various  density  and 
current  profiles  are  used  to  determine  the  dependence  of 
these  quantities  on  the  profile  shapes.  Our  results  demon¬ 
strate  the  importance  of  using  accurate  experimental  den¬ 
sity  and  current  profiles  with  radiation-scaling  analyses. 


I.  INTRODUCTION 

The  line  mass  m0  and  implosion  energy  Ew  inferred 
from  the  plasma  radiation  source  (PRS)  implosion  time 
and  load  current  depend  on  the  shapes  of  gas-puff  radial- 
density  and  current-history  profiles.  For  long  implosion 
times  and  high  photon  energy,  such  as  considered  for  the 
Decade  Quad  PRS  driver  [1],  computed  K-shell  yields  can 
depend  sensitively  on  these  quantities  [2].  Here,  snow¬ 
plow  calculations  with  various  density  and  current  pro¬ 
files  are  used  to  determine  the  dependence  of  m0  and  Ew 
on  their  shapes.  As  part  of  this  analysis,  a  ballistic-trans¬ 
port  model  is  developed  to  determine  realistic  density 
profiles  from  gas-puff-nozzle  parameters. 

Analytic  forms  relating  m0  to  the  density-  and  current- 
profile  shapes  are  found  when  a  mass-averaged  radius 
replaces  the  nominal  outer  radius  of  the  density  profile.  It 
is  shown  that  small  experimental  deviations  from  the 
commonly-used  linearly-rising  current  can  produce  fac- 
tor-of-two  errors  in  calculated  dynamic  mass,  and  that 
diffuse  density  profiles  (without  well-defined  current-ini¬ 
tiation  radii)  can  have  surprisingly-large  uncertainties  in 
implosion  energy.  Such  results  demonstrate  the  impor¬ 
tance  of  using  accurate  experimental  density  and  current 
profiles  with  radiation-scaling  analyses  [2,  3]. 


II.  CALCULATIONS  AND  DISCUSSION 

Radial  snowplow  implosions  of  a  PRS  with  initial  radial 
line-mass  distribution  m(r)  g/cm  driven  by  an  applied 
current  I(t)  A  are  described  by 
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where  R(t)  is  the  snowplow  radius, 
to 
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Scaling  Eq.  (1)  leads 
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where  tjmp  is  the  implosion  time,  Rc  is  a  characteristic  ra¬ 
dius,  and  I0  is  the  maximum  current  during  implosion. 
The  quantity  C,  varies  with  the  shapes  of  the  density  and 
current  profiles,  but  is  insensitive  to  the  value  of  stagna¬ 
tion  radius  Rf  =  R(timp).  Figure  1  shows  these  depen¬ 
dences  in  DM2  long-implosion-time,  argon-gas-puff  ex¬ 
periments  using  annular  and  solid-fill  nozzles  [4], 
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Figure  1.  Variation  of  implosion-time  factor  with 
implosion  time  for  three  DM2  nozzle  geometries. 
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In  Fig.  1,  m0  is  taken  to  be  proportional  to  the  product 
of  nozzle  pressure  p  and  throat  diameter  Dt  for  each  shot, 
while  the  mean  nozzle  diameter  D0  determines  R^  .  These 
data,  and  similar  data  from  experiments  carried  out  on 
HAWK  [5],  show  that  calculating  load  mass  from  the 
implosion  time  depends  on  the  density  distribution  and, 
through  its  variation  with  implosion  time,  current  shape. 

For  long-implosion-time  experiments  such  as  those  on 
DM2  and  HAWK  (where  a  weak  current  dip  at  stagnation 
is  observed),  a  linear  current  rise  to  a  flat  top  can  be  used 
to  characterize  current  profiles.  Designating  tc  as  the  time 
when  the  current  changes  from  rising  to  constant,  tc/timp  is 
in  the  range  0,65-0.85  for  DM2,  and  in  the  range  0.2-0. 6 
for  Hawk  operating  with  a  plasma  opening  switch  (POS). 

This  current  profile  is  first  used  with  a  variable-width, 
uniform-density  annulus  (UFA)  defined  by  n  =  constant 
for  Ri  <  R  <  Ro,  and  zero  otherwise.  Results  of  these 
snowplow  computations  can  be  summarized  by  two  mass 
factors  (~1/Ct2)  defined  by 
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is  a  mass-averaged  radius,  and  the  right-hand  side  of 
Eq.  (5)  is  for  the  UFA  density  distribution.  The 
variations  of  Cm  and  Cm'  with  annulus  width  A  =  R0  -  Rt 
are  shown  in  Fig.  2  for  a  linear  current  rise  (tc  =  timp).  The 
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Figure  2.  Variation  of  the  mass  factors  with  UFA- 
model  annulus  width  for  a  linear  current  rise. 

relative  constancy  of  Cm'  suggests  that  the  factor-of-two 
variation  in  m0  associated  with  density-profile  shape  can 
most-easily  be  accounted  for  by  using  Eq.  (4),  and  this 
form  is  displayed  in  what  follows.  Knowledge  of  the 
density  profile  is  still  required  to  determine  m0  from  <R>. 
This  procedure  was  used  in  HAWK  neon  gas-puff 
experiments  [5]  to  achieve  agreement  with  interferometer 
measurements  [6], 

Figure  3  plots  the  variation  of  Cm'  with  current-profile 
shape  for  three  values  of  A/Ro .  The  figure  shows  that  m0 
depends  sensitively  on  the  current  shape,  so  that  even  the 
modest  DM2  flat-top  current  can  increase  the  mass  in¬ 
ferred  from  a  linearly-rising  current  by  a  factor  of  2.  Cm' 
remains  nearly  constant  with  density-profile  shape.  The 
mass-factor  in  the  figure  is  well  fit  by 


Figure  3.  Variation  of  Cm'  with  current  shape  for 
various  UFA  density-profile  widths. 
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A  ballistic-transport  (BT)  model  can  be  used  to  deter¬ 
mine  more-accurate  density  profiles  from  gas-puff-nozzle 
parameters  [6].  For  this  model, 


n(r,z )  = 
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Texp 
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In  Eq.  (7),  N  is  the  line  density,  0d  is  the  divergence  angle 
of  gas  escaping  from  the  nozzle,  z  is  the  distance  from  the 
nozzle,  Rn  is  the  nozzle  radius,  0t  is  the  nozzle  tilt  angle, 
and  I0  is  the  Bessel  function.  Experimental  density  con¬ 
tours  measured  with  interferometry  at  various  distances 
from  the  nozzle  compare  favorably  with  those  predicted 
by  the  BT  model  [6]  using  single  values  for  the  three  noz¬ 
zle  parameters.  Partial  results  are  shown  in  Fig.  4. 

For  the  implosion  calculations  of  interest  here,  BT  den¬ 
sity  profiles  can  be  specified  in  terms  of  a  single  parame¬ 
ter  5/Ro  as  follows. 
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For  8/Ro  «  1,  the  radial  density  distribution  is  a  thin 
Gaussian  annulus  centered  about  R0 .  For  8/Ro  >  1,  i.e.  for 
z  >  R?s/(0d  +  0,),  the  profile  approaches  a  Gaussian  of  ex¬ 
panding  width  8  about  the  axis  of  symmetry.  Note  that 
for  BT  distributions,  R0  represents  a  radius  inside  the  dis¬ 
tribution  and  does  not  characterize  the  outside  edge, 
which  is  not  well-defined  for  such  diffuse  distributions. 
The  mass-averaged  radius  from  Eq.  (5)  is  well  fit  by 
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Figure  4.  Comparison  of  BT  model  (dashed)  and 
measured  (solid)  arial-density  distributions  at  three 
axial  locations. 


6 


4 

d 

2 


0 

0.0  0.5  1.0 

^c^imp 

Figure  5.  Variation  of  BT-model  mass-factor  with 
current  shape  for  various  8/Ro  values. 

similar  to  those  of  UFA  profiles,  though  the  variations  of 
m0  with  density  profile  through  Eq.  (4)  is  different  be¬ 
cause  of  the  different  dependence  of  <R>  on  radial  shape. 
For  BT  profiles,  <R>  >  R0  because  of  mass  contributions 
in  the  tail  at  larger  radius.  For  UFA  profiles,  <R>  <  R0 
since  all  mass  is  at  smaller  radius. 

The  S/Ro  =  0  curve  in  Fig.  5  is  the  same  as  that  of  the 
UFA  model  with  A/R0  =  0.  For  large  values  of  5/R0,  the 
shape  of  the  density  profile  remains  Gaussian  and  the  Cm' 
values  are  the  same  as  those  of  the  3.2  curve.  However, 
the  corresponding  m0  values  continue  to  change  because 
<R>  scales  with  5  as  in  Eq.  (9).  The  curves  of  Fig.  5  are 
well  fit  by  the  form  of  Eq.  (6)  with  12  +1.3A/Ro  replaced 
by  11.8  +3.5(8/Ro),  10  replaced  by  9.6  +1.9(8/Ro),  and  8 
replaced  by  7.5  +8/R0. 

For  snowplow  implosions,  the  implosion  energy  Ew 
exceeds  the  slug-model  kinetic  energy  and  can  be  deter¬ 
mined  from  the  electrical  characteristics  of  the  PRS  load. 
For  l(t)  in  MA  and  load  inductance  L(t)  in  nH/cm 

If  the  initial  radius  of  current  flow  Rj  is  well  defined,  Ew 
can  be  written  [2] 


■  (9) 
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The  mass  factor  Cm'  (8/Ro ,  te/timp)  can  be  determined  for 
these  more-realistic  density  profiles.  Related  phenomena, 
such  as  zippering,  can  then  be  determined  as  functions  of 
of  the  nozzle  parameters  Rn  ,  0d ,  and  0t . 

Figure  5  plots  the  variation  of  Cm'  with  current-profile 
shape  for  various  values  of  8/Ro.  BT  density-profile 
mass-factor  variations  with  current  shape  are  seen  to  be 


where  y  depends  on  the  density  and  current  shapes.  When 
y  =  1,  Eqs.  (10)  and  (11)  yield  I2AL/2  for  implosion  at 
constant  current  (tc  =  0).  Figure  6  shows  the  variation  in  y 
with  current  shape  for  UFA-model  density  distributions. 
The  implosion  energy  varies  by  less  than  10%  with  den¬ 
sity  distribution  and  can  be  determined  with  accuracy  for 
the  experimental  current  profile  [2].  Values  of  y  for  short- 
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Figure  6.  Variation  of  energy  factor  with  current 
shape  for  various  annular  width  in  the  UFA  model. 


implosion-time  current  profiles  with  strong  dips  near 
stagnation  will  be  smaller  by  about  0.1. 

For  diffuse  profiles  where  Rj  may  not  be  well-defined,  y 
may  be  replaced  by  ew  =  Ew(kJ/cm)/I02(MA),  computed 
from  Eq.  (10)  with  the  assumption  that  the  current  moves 
with  the  snowplow.  Results  for  BT  density  profiles  are 
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Figure  7.  Variation  of  implosion  energy  with  current 
shape  for  various  values  of  8/Ro  and  Rf  =  <R>/10. 

shown  in  Fig.  7  for  the  case  Rf  =  <R>/10.  Note  that  the 
compression  ratio  varies  with  5/R0  since  <R>  varies  with 
the  profile  shape.  For  UFA  profiles,  ew  increases  with  A 
because  compression  increases  like  10Ro/<R>,  and  <R> 
decreases  with  A.  For  the  diffuse  BT  profiles,  ew  in¬ 
creases  with  8  because  compression  increases  like 
10R;/<R>,  and  Rj  ,  far  out  on  the  tail  of  the  distribution, 
increases  with  8  faster  than  <R>  does.  The  BT  implosion 
energy  increases  at  small  tc/timp  because  the  current  is  es¬ 


tablished  at  early  times  when  the  snowplow  is  still  at  large 
radius.  Because  of  the  low-density  wings,  the  diffuse- 
profile  implosion  energy  depends  more  sensitively  on 
details  of  both  the  density  and  current  distributions.  This 
can  lead  to  large  uncertainties  in  the  calculation  of  Ew  if 
the  experimental  profiles  are  not  taken  into  account. 


III.  SUMMARY  OF  RESULTS 

The  dependence  of  dynamic  mass  mo  and  implosion 
energy  Ew  on  the  shapes  of  the  radial-density  and  current- 
history  profiles  has  been  systematically  quantified.  In 
addition  to  a  simple,  uniform-fill,  fat-annulus  density  pro¬ 
file,  a  general ly-useful,  ballistic-transport  .model  has  been 
developed  for  this  analysis  and  benchmarked  against  in¬ 
terferometer  measurements  to  determine  more-realistic 
density  profiles  from  gas-puff-nozzle  parameters.  Ana¬ 
lytic  forms  relating  m0  to  the  density-  and  current-profile 
shapes  are  found  when  a  mass-averaged  radius  replaces 
the  outer  radius  of  the  density  distribution.  Small  experi¬ 
mental  deviations  from  the  commonly-used  linearly-rising 
current  have  been  shown  to  produce  factor-of-two  errors 
in  dynamic  mass.  Also,  diffuse  density  profiles  without 
well-defined  current-initiation  radii  can  lead  to  surpris- 
ingly-large  errors  in  implosion  energy.  As  computed  K- 
shell  yields  can  depend  sensitively  on  m0  and  Ew  ,  our 
results  demonstrate  the  importance  of  using  accurate  ex¬ 
perimental  density  and  current  profiles  with  radiation¬ 
scaling  analyses. 
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